Pb-graphene-Pb Josephson junctions: characterization in B field. 
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We fabricate superconductor- graphene-superconductor Josephson junctions with superconducting 
regions made of lead (Pb) . The critical current through grapehene may be modulated by external 
magnetic field; the resulting Fraunhofer interference pattern shows several periods of oscillations, 
indicating that the junction is uniform. Deviations from the perfect Fraunhofer pattern are observed, 
and their cause is explained by a simulation that takes into account the sample design. 
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The properties of superconductor-graphene- 
superconductor (SGS) junctions have attracted 
significant attention [T]-[5]. Unlike the conventional 
superconductor-normal metal-superconductor (SNS) 
junctions, devices made with graphene allow for high 
tunability with the gate voltage. We have recently 
reported on SGS junctions which use lead (Pb) as the 
superconducting material [5 . Lead has a relatively high 
critical temperature of 7.2 K, compared to 1.2 K in the 
case of the commonly used aluminum. Indeed, we have 
observed supercurrent through graphene at tempera- 
tures as high as 2 K. In this paper, we characterize the 
properties of these junctions by applying magnetic field. 

We fabricate the superconducting contacts to graphene 
from a Palladium/Lead (Pd/Pb) bilayer. First, we de- 
posit a 2 nm layer of palladium, which creates transpar- 
ent contacts to graphene [6^ 7 ; lOOnm layer of lead is de- 
posited in situ on top. The lateral width of the contacts 
is typically 500 nm. In this work we present the results 
measured on a junction about 20/im wide and 400 nm 
long. (Commonly, length of the junction is defined as 
the distance between the superconducting contacts, and 
width as the distance of the normal metal along the su- 
perconducting contacts.) In order to create such a wide 
junction, the leads are bent in two places to fit on a 
moderately-sized graphene flake (Fig. la). We show that 
this particular sample design has certain nontrivial con- 
sequences. 

We measured the sample electronic properties using a 
pseudo 4-probe setup (Fig. la). The junction is biased 
by current / which contains a small AC component, and 
the AC voltage across the junction is measured using a 
lock-in amplifier. The carrier density in graphene can 
be tuned by the back-gate voltage Vgate but for the re- 
sults presented the gate voltage is set at zero. Finally, 
a perpendicular magnetic field can be applied using two 
methods. Conventionally, a field B^xt can be created by 
an external solenoid magnet. Alternatively, we send a 
current II along one of the superconducting leads (Fig. 
la), inducing a field which we label as Bl. The advan- 
tage of the second method is that the required small fields 
can be easily obtained and rapidly changed. In this sam- 
ple we have calibrated Bl to be equal to 0.95 ^/^ (see 



details below). 

The Pd/Pb electrodes become superconducting at a 
temperature of 7 K, and the SGS junctions begin to 
exhibit enhanced zero-bias conductance at temperatures 
of ^ 5 K. Below ^ 2 K, a fully formed supercurrent 
branch is clearly observed [5 . Figure Ic demonstrates 
the differential conductance R = dV/dl versus bias cur- 
rent / (vertical axis) and magnetic field Bl (horizontal 
axis) measured at several temperatures. The dark areas 
of the maps in Figure Ic correspond to the regions of sup- 
pressed resistance. The regions are bound by a critical 
current / = /c, above which the junction becomes nor- 
mal. The value of Ic increases as temperature is lowered 
and saturates around Ic ^ 0.5/iA at zero magnetic field 
(see the lowest map in Fig. Ic). When Bl is applied, 
Ic oscillates in a way closely resembling the Fraunhoffer 
diffraction pattern [10 . Several oscillations of Ic can be 
observed at the lowest temperature; this indicates that 
the junction is uniform. 

We next apply an external magnetic field B^xt , which is 
found to shift the modulation pattern of Figure Ic in the 
horizontal direction (Fig. 2). The shift is linear in B^xt'- 
indeed, at the center of the pattern the external field and 
the one induced by II cancel each other. The observed 
rate of shift allows us to fix the conversion Bl = 0.95 ^/l 
mentioned earlier. This factor is also consistent with our 
order of magniture estimates. Furthermore, the shift of 
the pattern by <l>o in an external magnetic field of 0.36 
mT allows us to extract the effective area of 5.6/im^. 
While this area is smaller than the Sfim? expected from 
the designed sample dimensions oi W = 20/im by L = 
OA/j^m, it is quite likely that the length L between the 
leads is reduced in the process of lithography, or that 
the magnetic field is modified due to the presence of the 
superconducting leads. 

When magnetic field B^^t of the order of tens of mT is 
applied to the sample, the observed pattern becomes dis- 
torted even after the field is returned back to zero (Fig. 
3). It is clear that the resulting pattern at B^^t = 
(Fig. 3a) is very different from the original one (Fig. 2c). 
We can partially recover the original pattern by setting 
Bext ^ 3.4mT (Fig. 3b). When comparing the resulting 
pattern to the original one (Fig. 2c), we notice that the 
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FIG. 1: b) Scanning electron micrograph of a graphene- 
based SNS Josephson Junction similar but smaller than the 
one presented in this paper. (Junctions used for measurement 
were not imaged in order to preserve the quality of graphene.) 
The dark triangular area is a single layer graphene flake, and 
the metal contacts are made from lead (Pd) with a thin con- 
tact layer of palladium (Pd). a) Schematic of the measure- 
ment setup. The metal leads form a "|J" shape in order to 
increase their length. Bias current / is sent through the junc- 
tion with a small AC component. The resulting AC compo- 
nent of the voltage across the junction is measured using a 
lock-in amplifier allowing one to record the differential resis- 
tance R = dV/dl. An external magnetic field Bext is applied 
by a superconducting solenoid. In addition, a magnetic field 
5l is created by sending a current II along one of the Pb 
leads of the junction. Sweeping the current II allows to ap- 
ply a very small magnetic field Bl- c) Differential resistance 
dV/dl maps measured ^;s. bias current / and magnetic field- 
inducing current II- Regions of vanishing R appear dark. 
Each panel corresponds to measurement at a different tem- 
perature. Enhanced zero-bias conductance develops around 
^ 4:K for small fields Bl- With lower temperatures more and 
more critical current modulations appear and the Fraunhofer 
interference pattern is observed. At the base temperature of 
1.3K, critical current is seen at fields beyond 5mT (only fields 
up to lAmT are shown). Observing many oscillations suggest 
that the junction is highly uniform. 
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FIG. 2: R(I,Il) maps (like those in Fig. Ic) measured at 
different values of the external magnetic field (a): Bext = 
-1.3 mT, (b): Bext = -0.36 mT, (c): Bext = 0, (d): Bext = 
0.36 mT, (e): Bext = 1.3 mT. Application of Bext shifts the 
modulation pattern, so that at its center Bext and Bl cancel 
each other (b-d). Since the cancellation is not perfect, the 
pattern gets distorted (a), compared to the pattern at zero 
external field (c). An opposite orientation of the external 
magnetic field (e) results in mirror reversal of the distortions. 
T= 1.3 K. 

critical current is slightly suppressed and the side-lobes 
have somewhat random heights. We attribute these dis- 
tortions and the shift from zero field to the trapping of 
magnetic flux in the superconducting film ^J. Indeed, 
the undistorted pattern seen in Figure 2c can be restored 
following the thermal cycling to ~ lOK, beyond the crit- 
ical temperature of lead. 

Interestingly, at fields less than those causing trapped 
flux of Figure 3, distortions of a different nature are in- 
troduced to the pattern (Fig. 2a). The first difference 
is that, returning to B^xt = restores the original pat- 
tern without any hysteresis. Second, the pattern demon- 
strates perfect symmetry under simultaneous reversal of 
both Bext and J^, compare Figures 2 (a) and (e). We 
associate this behavior with the fact that the phase dif- 
ference is not linear along the length of the leads. Indeed, 
the field Bl may not be entirely uniform, so that it is 
not perfectly compensated by B^xt- Most likely, the de- 
viations of Bl from uniformity are caused by the bends 
in the leads (schematic in Fig. la), at which points the 
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FIG. 3: R{I, II) measurements taken after the perpendicular 
magnetic field was ramped beyond several tens of mT. Panel 
(a) shows the measurement done when Bext was returned back 
to zero. Clearly, the pattern is now vastly distorted. Panel (b) 
is taken at Bext — 3.4mT. At this field, the original critical 
current modulation pattern (Fig. 2c) is partially restored. 
However, the central lobe shows a suppressed critical current, 
and the side-lobes form a distorted pattern. These permanent 
distortions are attributed to trapped flux in the Pb leads. 
Heating of the sample beyond the Tc of Pb is required in 
order to restore the symmetric patterns seen in Figure 2c. 



phase difference experiences discontinuous jumps propor- 
tional to The situation is very similar to the junctions 
with an artificial phase discontinuity controlled by an ex- 
ternal current [121 [13]. Indeed, some of the features we 
observe in Figures 2(a) and (e), e.g. the strengthening of 
the side lobe at the expense of the central lobe, resemble 
those found in Refs. [T21 IT3]. 

To describe the distortions found in Figure 2a,e, we 
consider a semi-realistic model of the sample. We assume 
that the leads extend from x = —W/2 to +1^/2; the 
position-dependent phase difference ^l(^) (induced by 
Bl{x)) is taken to be piece- wise linear in x, with a slope 
proportional to II- Two identical discontinuous jumps of 
are placed at -W/10 and ^W/10. The strength 
of the discontinuities is taken to be also proportional to 
II' These points are close to the actual locations of the 
bends in the leads, but we checked that the main features 
of the simulation do not crucially depend on the details 
(i.e. the position of discontinuities or symmetry of their 
placement). 

We also include the effect of the external field, pre- 
suming it induces uniform phase difference ^ext{x) along 
the length of the leads. In our simulations, the current- 
phase relation is assumed to be sinusoidal. Although 
deviations from a sinusoidal relation have been recently 
observed in SGS junctions [14 , the approximation should 
be adequate in our case, due to the relatively large length 
between the leads {L = 400 nm) and the relatively high 
temperature > IK. 

The simulated patterns of the critical current Ic vs. 
Il and Bext are shown in Figure 4. The main features 



observed in Figure 2 are qualitatively reproduced, such 
as: 1) the overall shift of the Ic{Il) pattern in B^xu 2) 
the growing distortion of the Ic{Il) pattern in Bext] 3) 
the growing strength of the side lobe on the high current 
side of the pattern at the expense of the central lobe; 

4) the difference in width between that side lobe and the 
side lobes on the other side. Moreover, even the Bext = 
curve, while similar to the perfect Fraunhofer pattern / ex 
smiirlL/ 1^^) / Ili bears noticeable differences. Namely, 
some of the side lobes are suppressed almost to zero, 
while further lobes at higher 1^ regain strength. This 
type of behavior is indeed observed in experiment (Fig. 

5) . Note the region of suppressed critical current at II ^ 
±2.3m74 in Figure 2c, and its reappearance at higher 
Il - ±3.5mA. 
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FIG. 4: Simulated critical current vs. Il at several values of 
Bext- The bottom curve is taken at Bext — 0, and in each con- 
secutive curvethe external flux grows in units of flux quantum 
$0- The phase difference between the two leads induced by 
Il is assumed to grow linearly along the length of the leads, 
proportionally to /l, and to jump discontinuously at two lo- 
cations by an amount also proportional to II • This particular 
functional form it chosen to approximate the realistic shape 
of the sample, where the leads turn 90° in two places; however 
the major features appear insensitive to the exact locations 
of the discontinuities. The horizontal axis is labeled in units 
of total flux ^L induced by II excluding the discontinuities. 
The additional phase jump at each discontinuity is equal to 
O.Itt^l/^o. 

In conclusion, we demonstrate magnetic field-induced 
quasi-periodic modulation of critical current in Pb- 
graphene-Pb structures, which indicates their spatial uni- 
formity. The magnetic field can be applied by running a 
current through one of the superconducting leads within 
the same structure, resulting in a simple, yet efficient 
method to in situ control the critical current. The de- 
pendence of the critical current on thus applied magnetic 
field deviates from the perfect Fraunhofer interference 
pattern. The difference is attributed to the presence of 
bends in the superconducting leads; a simple simulation 
supports this explanation. 
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FIG. 5: Differential resistance dV/dl map measured ^;s. bias current / and current II inducing magnetic field. This data are 
similar to Fig. 2c, but is extended up to II — QmA. Note the suppression of select side-lobes (see II ~ ±2. 5mA and ±3. 5mA). 
This suppression of side- lobes qualitatively resembles the simulation results (Fig. 4, Bext = 0, $l/^o ~ ±3 and ±5). The tilt 
of the pattern is an artifact of II flowing through the normal part of the sample, thus creating a voltage drop that shifts the 
zero of /. 
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